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^ ' Abstract 

a.: 



With its unprecedented light-collecting area for night-sky observations, 
the Cherenkov Telescope Array (CTA) holds great potential for also optical 
stellar astronomy, in particular as a multi-element intensity interferometer for 
realizing imaging with sub-milliarcsecond angular resolution. Such an order- 
of-magnitude increase of the spatial resolution achieved in optical astronomy 
will reveal the surfaces of rotationally flattened stars with structures in their 
^ . circumstellar disks and winds, or the gas flows between close binaries. Image 

reconstruction is feasible from the second-order coherence of light, measured 
as the temporal correlations of arrival times between photons recorded in 
different telescopes. This technique (once pioneered by Hanbury Brown and 
C\l ■ Twiss) connects telescopes only with electronic signals and is practically in- 

sensitive to atmospheric turbulence and to imperfections in telescope optics. 
Detector and telescope requirements are very similar to those for imaging 
air Cherenkov observatories, the main difference being the signal processing 
(calculating cross correlations between single camera pixels in pairs of tele- 
scopes). Observations of brighter stars are not limited by sky brightness, 
permitting efficient CTA use during also bright-Moon periods. While other 
concepts have been proposed to realize kilometer-scale optical interferometers 
of conventional amplitude (phase-) type, both in space and on the ground, 
their complexity places them much further into the future than CTA, which 
thus could become the first kilometer-scale optical imager in astronomy. 
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1. Resolution frontiers in astronomy 

Many efforts in optical astronomy aim at improving the spatial resolution 
in order to obtain ever sharper views of our Universe. Projects include the 
construction of extremely large telescopes utilizing adaptive optics, or plac- 
ing instruments in space. The highest resolution is currently obtained from 
amplitude (phase-) interferometers which combine light from telescopes sep- 
arated by baselines up to a few hundred meters. Since effects of atmospheric 
turbulence are less severe at longer wavelengths, such instruments are prefer- 
entially operated in the near infrared. Tantalizing results from such facilities 
show how stellar disks start to become resolved, beginning to reveal stars as 
a vast diversity of individual objects, although so far feasible only for a small 
number of the largest ones which extend for tens of milliarcseconds (mas). 
More typical bright stars have diameters of only a few mas, requiring inter- 
ferometry over many hundreds of meters or some kilometer to enable surface 
imaging. Using a simple A/r criterion for the required optical baseline, a 
resolution of 1 milliarcsecond at A 500 nm requires around 100 meters, while 
1 km enables 100 //as. 

Since we currently are at the threshold of starting to resolve stars as ex- 
tended objects, a great step forward will be enabled by improving the resolu- 
tion by just another order of magnitude. However, since ordinary amplitude 
interferometers require precisions in both their optics and in the atmosphere 
above to within a small fraction of an optical wavelength, atmospheric tur- 
bulence constrains their operation to baselines not much longer than some 
100 m, especially at shorter visual wavelengths. 

The scientific promise of very long baseline optical interferometry for 



imaging stellar surfaces has been realized by several [63|, [99|, and concepts 
to circumvent atmospheric turbulence include proposals for large amplitude 
interferometer arrays in space: Stellar Imager [rj| and the Luciola hypertele- 



scope [65J , or possi bly placed at extreme terrestrial locations such as Dome C 



in Antarctica [116J. However, the complexity and likely cost of these projects 
make the timescales for their realization somewhat uncertain, prompting 
searches for alternative approaches. Although not a complete replacement 
for the many capabilities of large space-based interferometers, comparable 
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science can begin to be realized very much sooner, and with much less ef- 
fort, by ground-based intensity interferometry, utilizing large arrays of air 
Cherenkov telescopes. 

1.1. Intensity interferometry 

Intensity interferometry was pioneered by Robert Hanbury Brown and 



Richard Q. Twiss 35], for the original purpose of measuring stellar sizes, 
and a dedicated instrument was built at Narrabri, Australia. It measures 
temporal correlations of arrival times between photons recorded in different 
telescopes to observe the second-order coherence of light (i.e., that of in- 
tensity, not of amplitude or phase). The name intensity interferometer is 
sort of a misnomer: actually nothing is interfering in the instrument. This 
name was chosen for its analogy to the ordinary amplitude interferometer, 
which at that time had similar scientific aims in measuring source diameters. 
Two separate telescopes simultaneously measure the random and very rapid 
intrinsic fluctuations in the light from some particular star. When the tele- 
scopes are placed sufficiently close to one another, the fluctuations measured 
in the two telescopes are correlated, but when moving them apart, the fluc- 
tuations gradually become independent and decorrelated. How rapidly this 
occurs for increasing telescope separations gives a measure of the spatial co- 
herence of starlight, and thus the spatial properties of the star. The signal is 
a measure of the second-order spatial coherence, the square of that visibility 
which would be observed in any classical amplitude interferometer. Spatial 
baselines for obtaining any given resolution are thus the same as would be 
required in ordinary interferometry. 

The great observational advantage of intensity interferometry (compared 
to amplitude interferometry) is that it is practically insensitive to either 
atmospheric turbulence or to telescope optical imperfections, enabling very 
long baselines as well as observing at short optical wavelengths, even through 
large airmasses far away from zenith. Telescopes are connected only with 
electronic signals (rather than optically), from which it follows that the noise 
budget relates to the relatively long electronic timescales (nanoseconds, and 
light-travel distances of centimeters or meters) rather than those of the light 
wave itself (femtoseconds and nanometers). A realistic time resolution of 
perhaps 10 nanoseconds corresponds to 3 m light-travel distance, and the 
control of atmospheric path-lengths and telescope imperfections then needs 
only to correspond to some reasonable fraction of those 3 meters. 
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The measured second-order coherence provides the square of the ordinary 
visibility and always remains positive (save for measurement noise), only di- 
minishing in magnitude when smeared over time intervals longer than the 
optical coherence time of starlight (due to finite time resolution in the elec- 
tronics or imprecise telescope placements along the wavefront). However, for 
realistic time resolutions (much longer than an optical coherence time of per- 
haps ~10~ 14 s), the magnitude of any measured signal is tiny, requiring very 
precise photon statistics for its reliable determination. Large photon fluxes 
(and thus large telescopes) are therefore required; already the flux collectors 
used in the original intensity interferometer at Narrabri were larger than any 
other optical telescope at that time. 

Details of the original intensity interferometer at Narrabri and its observ- 
ing program (mainly measuring angular sizes of hot stars) were documented 



by Hanbury Brown et al. [39|, |40j, including retrospective over views |35l 



38| . The principles are also explained in various textbooks [64], llOOL 1108 . 

The original intensity interferometer at Narrabri had two reflecting tele- 
scopes of 6.5 m diameter, formed by mosaics of numerous hexagonal mirrors, 
providing star images of 12 arcmin diameter. Following the completion of 
that program, the design for a second-generation intensity interferometer 



was worked out [16|, |36|, l38[. This larger facility was envisioned to have 12-m 
diameter telescopes, movable over 2 km, however it was never realized. The 
same physical principles of measuring intensity correlations have since been 
actively utilized in high-energy particle physics (where also other bosons, 
i.e., particles with an integer number of their quantum spin, have a ten- 
dency to bunch together in a similar way as photons, while electrons and 
other fermions show the opposite behavior). In astronomy, however, inten- 
sity interferometry has not undergone further development, largely due to 
its demanding requirements for large optical flux collectors, spread over long 
baselines, and equipped with fast detectors and high-speed electronics. 



1.2. Air Cherenkov telescopes 

The parameters of air Cherenkov telescopes are remarkably similar to the 
requirements for intensity interferometry. In the Narrabri interferometer, 
movable telescopes were used to maintain a fixed baseline while tracking 
a source across the sky. Nowadays, electronic time delays can compensate 
for the different arrival times of a wavefront to different telescopes in fixed 
positions. 
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The most remarkable potential comes from the Cherenkov Telescope Ar- 
ray [H, [3| which foresees a total of 50-100 telescopes with differently sized 
apertures between about 5 and 25 meters, distributed over an area of 2- 
3 km 2 . Such a large array permits an enormous number of optical baseline 
pairs to be synthesized, enabling measurements of angular scales between 
milli- and microarcseconds. The potential of using such arrays for inten- 



sity interferometry has indeed been noticed by several authors [19|, [69|, 171 
Within the CTA project, a task group was set up to specify how to enable 
it for also such uses. If a baseline of 2 km could be utilized at A =350nm, 
resolutions would approach 30 //as, an unprecedented spatial resolution in 
optical astronomy. Such numbers are challenged only by radio interferome- 



ters operating between Earth and antennas in deep space [56j, or possibly by 



futuristic X-ray interferometers [78 



2. Principles of intensity interferometry 

In its simplest form, an intensity interferometer consists of two optical 
telescopes or light collectors, each with a photon detector feeding one chan- 
nel of a signal processor for temporally cross correlating the light-intensity 
signals from the two telescopes. The intensities measured at detectors 1 and 
2 are the respective values of the electric light- wave amplitude times its com- 
plex conjugate, averaged over some time interval corresponding to the signal 
bandwidth of the detectors and associated electronics: 



(1(f)) = (E(t)E*(t)} (1) 

where * marks complex conjugate and ( ) denotes averaging over time. 
The intensities measured in the two telescopes are cross correlated: 



(h(t)I 2 (t)) = (Ei(t)E{(t) ■ E 2 (t)E*(t)) (2) 

This expression can be expanded by dividing the complex field amplitudes 
into their real and imaginary parts. Here one must make an assumption that 
is fundamental to the operation of an intensity interferometer: the light must 
be chaotic, i.e., with a Gaussian ampli tude distribution; also called thermal- 



or maximum-entropy light [4|, |3l|, |73j, Il08| . To a good approximation this 



applies to all 'ordinary' light sources (but not necessarily to nonthermal ones 
such as lasers). Such light may well be quasi- monochromatic, as long as the 
light waves undergo random phase shifts, so that intensity fluctuations result 



5 




Figure 1: Principle of a multi-element stellar intensity interferometer. Several telescopes 
observe the same source, simultaneously recording its rapidly fluctuating optical light 
intensity I n {t). Cross correlations of the intensity fluctuations are measured between 
different pairs of telescopes: (h(t)I 2 (t)), (h(t)I 3 (t)), (h(t)h(t)), (I 2 (t)I 3 (t)), (I 2 (t)h(t)), 
etc. These yield a measure of the second-order spatial coherence of light, from which an 
image of the source can be deduced, with an angular resolution corresponding to the 
optical diffraction over the projected baseline distance between each pair of telescopes. 
Numerous telescopes enable a very large number of baselines to be synthesized, permitting 
a high-fidelity reconstruction of the source image. Telescopes distributed over also km-long 
baselines enable an angular resolution so far unprecedented in optical astronomy. 
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on timescales corresponding to the optical coherence time. For chaotic light, 
the real and imaginary parts of E\ and E 2 are Gaussian random variates, i.e., 
the values of E% and E 2 measured at different times can be treated as random 
variables obeying a normal distribution. Then the Gaussian moment theorem 
applies, which relates all higher-order correlations of Gaussian variates to 
products of their lower-order correlations (described in detail by Mandel & 
Wolf H)- 

It is then possible to show j64| that, for linearly polarized light: 



(J 1 (t)/ 2 (t)) = (/ 1 )(/ 2 )(l + | 7l2 | 2 ) (3) 

where 712 is the mutual coherence function of light between locations 1 
and 2, the quantity measured in ordinary amplitude interferometers. 
Defining the intensity fluctuations AJ as: 

Ah(t) = h(t) - (h) A/ 2 (t) = I 2 (t) - (I 2 ), 
one obtains: 

(A/ 1 (t)A/ 2 (t)) = (/ 1 )(/ 2 )| 7l2 | 2 , (4) 

since (A/) = 0. 

An intensity interferometer thus measures |7i 2 | 2 with a certain electronic 
time resolution. This quantity remains positive irrespective of atmospheric 
or optical disturbances although - since realistic time resolutions do not 
reach down to optical coherence times - it may get strongly diluted relative 
to the full value it would have had in the case of a hypothetical 'perfect' 
temporal resolution (shorter than the light- wave period). For realistic values 
of nanoseconds, this dilution typically amounts to several orders of magnitude 
and thus the directly measurable excess correlation becomes quite small. This 
is the reason why very precise photon statistics are required, implying large 
flux collectors. 



3. Optical aperture synthesis 

The original intensity interferometer at Narrabri used two telescopes, 
movable on railroad tracks, which could be positioned at different separa- 
tions r, to deduce angular sizes of stars from the observed function |7i 2 (r)| 2 , 
analogous to what can be measured with a two-element amplitude interfer- 
ometer. Systems with multiple telescopes and different baselines (Figure [1]) 
enable correspondingly more complete image reconstructions. Techniques for 
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interferometric imaging and aperture synthesis were first developed for ra- 



dio tele scop es [11 Ik Ill2j , but have since been elaborated also for the optical 
34j . 164 llOOj . Here we recall the basics: 

The separation vector between a pair of telescopes in a plane perpendic- 
ular to the line of observation, the (u,v)— plane, is ri — r2, so that for an 
optical wavelength A, r x — r 2 = (u\,v\). If the telescopes are not in such a 
plane, also a third coordinate enters: the time-delay w for the propagation 
of light along the line of sight to the source; ri — r 2 = (u\, v\, w). 

With the angular coordinate positions of the target (7, m), one can deduce 
the following expression for the correlation function Ti2 = (E(ri)E*(r2)): 

r M =//,(,, m) e~— (5) 

This equation represents the van Cittert-Zernike theorem, equating the 
quantity measured by an [amplitude] interferometer for a given baseline to 
a component of the Fourier transform of the surface intensity distribution of 
the source. This Fourier transform can be inverted: 

I v (l,m)= 1 1 ' V(u,v)e 2ni{ul+vm) dudv, (6) 



where V(u,v) equals the normalized value of j(u, v). Thus, by using 
multiple separations and orientations of interferometric pairs of telescopes, 
one can sample the (u,v)— plane and reconstruct the source image with a 
resolution equal to that of a telescope with a diameter of the longest baseline. 
This is the technique of aperture synthesis. 

In intensity interferometry, however, an additional complication enters 
in that the correlation function for the electric field, 712, is not measured 
directly, but only the square of its modulus, | ^yi2 1 2 - Since this does not 
preserve phase information, the direct inversion of the above equation is not 
possible. 

This limitation will be removed in intensity interferometry carried out 
with larger telescope arrays. For CTA, with some 50 or more elements, the 
possible number of baselines between telescope pairs becomes enormous; N 
telescopes can form N(N — 1)/2 baselines, reaching numbers in the thousands 
(even if possibly periodic telescope locations might make several of them re- 
dundant). Since such telescopes are fixed on the ground, the projected base- 
lines trace out curves in the (u,v)— plane, as a source moves across the sky. 
With proper signal handling, all successive measures of I712I 2 can be allocated 
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to their specific (u, v)— coordinates, producing a highly filled (u, v)— plane, 
with a superior coverage of projected orientations across the source image. 
As will be discussed below, such complete data coverage indeed enables re- 
construction of the phases of the Fourier components, and thus permits full 
two-dimensional image reconstructions (although the completeness of such 
coverage depends on how the source moves across the sky and thus on, e.g., 
whether it is located near the celestial equator or close to its poles). 

For large numbers of telescopes, another advantage of intensity inter- 
ferometry becomes obvious. Since telescopes connect only with electronic 
signals, there is in principle no loss of data when synthesizing any number 
of baselines between any pairs of telescopes: the digital signal from each 
telescope is merely copied electronically. By contrast, amplitude interferom- 
etry in the optical (as opposed to radio) requires optical beams of actual 
starlight between telescopes since the very high optical frequency (combined 
with rapid phase fluctuations in chaotic light) precludes its amplification 
with retained phase information. In order to obtain the many baselines 
needed for efficient aperture synthesis (such as realized in radio), starlight 
from each telescope is split and sent to beam combiners to interfere with the 
light from other telescopes, each combination with its own delay-line system. 
While such ambitious arrangements can be made for a moderate number 
of telescopes [l3j], the complexity (and the dilution of light between differ- 
ent baselines) rapidly increases if any greater number of telescopes would be 
engaged. 



4. The Cherenkov Telescope Array 

CTA is envisioned to have on the order of 50-100 telescopes with vari- 
ous apertures between about 5 and 25 meters, with currently favored con- 
figurations reaching an edge-to-edge distance of some 2 km. A number of 
candidate array layouts for the CTA were considered within its design study 
[3,0, 14, 44]: Figure [21 of which examples representing qualitatively different 



types of layouts are in Tabled] For interferometry, large telescope separations 
(long baselines) measure high-frequency Fourier components, corresponding 
to small structures on the target, while short baselines sample the low fre- 
quencies. For an Earth-bound interferometer (in a plane perpendicular to 
the line of observation) with a baseline B = (-£>North, -^East) the associated 
coordinates in the Fourier (u,v)— plane are (u, v ) = -J-(£>North, -BEast)- 
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Figure 2: Different telescope array layouts evaluated within the CTA design study, and 
also considered for their suitability to intensity interferometry. Each of the configurations 
labeled A through K is a subset of the all-encompassing hypothetical large array shown 
at bottom right. In this work, configurations B, D and I were selected as representative 
for three different classes of array geometry. 
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Figure 3: Coverage of the interfcrometric (u,v)— plane for three types of evaluated CTA 
layouts (top to bottom: B, D and I of Figure [2]). Left: (u,v)— plane coverages at one 
instant in time, for a star observed in the zenith. Upper right-hand squares expand the 
central 400x400m area. Right: (u,v)— plane coverages for a star moving from the zenith 
through 20 degrees to the west. The numerous telescopes enable a huge number of baseline 
pairs which largely fill the entire (it, v)— plane. 
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For stationary telescopes, the projected baselines, B p , will change while 
the target of observation moves across the sky, with each telescope pair trac- 
ing out an ellipse in the Fourier plane according to the following expression 
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where / is the latitude of the telescope array, and 5 and are the declina- 
tion and hour angles of the star. The w component corresponds to the time 
delay in the wavefront arrival time between the two telescopes (dependent 
on also the elevation difference of the telescopes, B up ). The extensive cov- 
erage of the (u,v)— plane that results from the Earth's rotation enables the 
synthesis in software of a very large telescope and - of course - is the very 
principle used in much of radio interferometry. 

Figure |3] illustrates these capabilities for three among the potential lay- 
outs considered for CTA, here taken as examples of qualitatively different 
telescope arrangements. One is a compact configuration; another a sparse 
and rather uniform one; and a third has telescopes of different sizes grouped 
with successively different spacings. The latter type of layout seems to lie 
close to those currently favored for CTA in general, and is also the most ca- 
pable one for interferometry. As seen in Figure |3j already short observations 
of just an hour or so, may cover much of the (u, v)— plane (and the coverage 
can be increased by observing at different wavelengths). 



Table 1: Properties of the three examined array layouts (B, D, and I in Figure^) from the 
CTA design study N is the number of telescopes, A is the light-collection area of each 
type of telescope, b is the number of unique baselines available, B m i n , B max indicate the 
range of baselines for observations in zenith. The corresponding range of angular diameters 
in milliarcseconds (1.22A/V) for observations at A 400 nm is indicated by 6 m in,dmax- 
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5. Signal-to-noise ratios in intensity interferometry 



For one pair of telescopes, the signal-to-noise ratio is given by 35|, Ill4j : 



(S/N) RMS = A-a-n- | 7 i 2 (r)| 2 • A/ 1/2 • (T/2) 1 / 2 (8) 

where A is the geometric mean of the areas (not diameters) of the two 
telescopes; a is the quantum efficiency of the optics plus detector system; n 
is the flux of the source in photons per unit optical bandwidth, per unit area, 
and per unit time; |7i2(r)| 2 is the second-order coherence of the source for 
the baseline vector r, with 712(1") being the mutual degree of coherence. A/ 
is the electronic bandwidth of the detector plus signal-handling system, and 
T is the integration time. 

Most of these parameters depend on the instrumentation, but n depends 
on the source itself, being a function of its radiation temperature. For a given 
number of photons detected per unit area and unit time, the signal-to-noise 
ratio is better for sources where those photons are squeezed into a narrower 
optical band. This property implies that (for a flat-spectrum source) the S/N 
is independent of the width of the optical passband, whether measuring only 
the limited light inside a narrow spectral feature or a much greater broad- 
band flux. Although perhaps somewhat counter-intuitive, the explanation is 
that realistic electronic resolutions of nanoseconds are very much slower than 
the temporal coherence time of broad-band light (perhaps 10~ 14 s). While 
narrowing the spectral passband does decrease the photon count rate, it also 
increases the temporal coherence by the same factor, canceling the effects of 
increased photon noise. This property was exploited already in the Narrabri 



interferometer [41| to identify the extended emission-line volume from the 
stellar wind around the Wolf-Rayet star 7 2 Vel. The same effect could also 
be exploited for increasing the signal-to-noise by observing the same source 
simultaneously in multiple spectral channels, a concept foreseen for the once 



proposed successor to the original Narrabri interferometer [16|, |36|, [38 . 



6. Simulated observations in intensity interferometry 

To obtain quantitative measures of what can be observed using realistic 
detectors on Cherenkov telescopes, a series of simulations were carried out. 
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6.1. Numerical simulations 

An intensity interferometer using two photon-counting detectors A and B 
and a digital correlator measures the squared modulus of the complex degree 
of coherence of the light: 

l2 (Aii A/2) (g) 



<w 2 > 



or, in a discrete form: 



^ = & < 10 > 

where iV^ and Nb are the number of photons detected in A and B re- 
spectively, Nab is the number of joint detections (i.e., the number of time 
intervals in which both detectors record a photon), and N is the number of 
sampled time intervals. Since a strict Monte-Carlo simulation would be com- 
putationally very demanding, a simplified procedure was used by generating 
random numbers Na, Nb and Nab, and inserting these into Eq. ([TO]) . These 
are Poisson-distributed random variables with mean values \x& = Pa • N, 
Hb = Pb ■ N and \lab = Pab • N. Here, Pa and Pb are the probabilities of 
detecting a photon in A and B respectively, within a small time interval At, 
and Pab is the probability of a joint detection within At. 

These probabilities can be written out in terms of variables depending 
only on the instrumentation and the target of study: 

Pa = a A {I A )At (11) 
P B = a B (I B )At (12) 
Pab = PaPb + a A a B {I A ){lB)\lAB\ 2 r c At (13) 

Here a denotes the quantum efficiency of the detectors, (J) is the mean 
light intensity, r c is the coherence time of the light (determined by the wave- 
length and optical passband) and jab is the degree of optical coherence 
(proportional to the Fourier transform of the target image, assuming tele- 
scope sizes to be small compared to the spatial structure in this transform). 
Such simulations were carried out for various telescope-array configurations 
and for various assumed sources. Here, examples are shown for a close binary 
star with components taken as uniform disks of diameters 200 and 150 /xas. 
Both the pristine original image and its pristine Fourier transform in the 
(u, v)— plane are shown in Figure |H Across the Fourier plane, the magnitude 
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Figure 4: Image of a close binary star with 200 and 150 //as diameter components, used to 
simulate observations, and the (logarithmized) magnitude of its Fourier transform. This 
noise-free pattern is what would be measured by a perfect interferometer of projected 
size 2000 x2000 m. Corresponding patterns in later figures cover only some part of this 
(it, v)— plane (due to finite extent of the telescope array on the ground), and become noisy 
for fainter sources and finite integration times. 



of various patterns varies greatly. To enhance the visibility of also fainter 
structures (and later to better see the effects of noise), the Fourier-plane 
figures use a logarithmic scaling and a shading to enhance the contrast (the 
exact numerical values of the measured correlations are not significant in this 
context). 

Also results from the simulated observations are mostly given as such 
Fourier-plane images rather than full image reconstructions. The simu- 
lated observations produced values at many different discrete locations in 
the (u,v)— plane, which were used in a linear interpolation to obtain the 
Fourier magnitude over a regular grid. This image format makes the effects 
of noise and changing telescope arrangements easier to interpret since it is 
independent of the performance of algorithms for image reconstruction or 
data analysis. As discussed below, optimal image reconstruction is a devel- 
oping research topic of its own. Even though reconstructed images do reflect 
the capability of the simulated telescope array, some reconstructions are still 
limited by the algorithms used. By contrast, the information recovered in 
the (u,v)— plane is independent of such algorithm performance. 

6.2. Limiting stellar magnitudes 

The question of how faint sources that can be usefully observed has been 



examined [28|, [7l|, with the conclusion that a conservative practical limit 
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for two-dimensional imaging with a large array of the CTA type is around 
niy=6. However, if only some one- dimensional measure would be sought 
(e.g., a stellar diameter or limb darkening), the data can be averaged over all 
position angles, and the limiting magnitude will become somewhat fainter. 
In any case, there are thousands of stars bright enough to be observable. 



7. Imaging with intensity interferometry 

An intensity interferometer directly measures only the absolute magni- 
tudes of the respective Fourier transform components of the source image 
that cover the (u, v)— plane, while the phases are not directly obtained. Such 
Fourier magnitudes can well be used by themselves to fit model parameters 
such as stellar diameters, stellar limb darkening, binary separations, circum- 
stellar disk thicknesses, etc., but two-dimensional images cannot be directly 
computed from the van Cittert-Zernicke theorem, Eq.([6]). However, a multi- 
component interferometer offers numerous baselines, and gives an extensive 
coverage of the (u, v)— plane, and it is already intuitively clear that the infor- 
mation contained there must place rather stringent constraints on the source 
image. 

7.1. Phase reconstruction 

A number of techniques have been developed for recovering the phase of 
a complex function when only its magnitude is known. Methods specifically 



intended for intensity interferometry have been worked out for one [49 



or 



two dimensions [50j. Once a sufficient coverage of the Fourier plane is avail- 
able, and phase recovery has been performed, image reconstruction becomes 
straightforward. 

7.2. Image reconstruction 

Various mathematical methods (of different numerical sophistication and 
sensitivity to various types of noise) can be applied for the reconstruction 
of images, and the development of optimum algorithms is a research topic 
of its own like, perhaps, was the case in early radio interferometry, before 
today's standard procedures (such as CLEAN) were developed. Nunez et 
al. (sESS 

applied Cauchy-Riemann based phase recovery to reconstruct 



images from simulated observations of oblate rotating stars, binary stars, and 
stars with brighter or darker regions, demonstrating that also rather complex 



16 



Reconstructed 




-1.0-0.5 0.0 0.5 1.0 1.5 
Milliarcseconds 



-1.0-0.5 0.0 0.5 1.0 
Milliarcseconds 



1.0-0.5 0.0 0.5 1.0 
Milliarcseconds 



|0.4 t 



Figure 5: Reconstructed images of binary stars (with varying diameter of the secondary) 
from simulated CTA observations. These simulations were for the array layout B (Fig- 
ure [2]), for sources assumed to have visual magnitude my = 3, and effective temperature 
T e ff = 7000 K. The assumed pristine images are shown below while the corresponding 
(u, v)— plane coverage is in Figure [H 




u-A (m) u-A (rrt) u-A (m) 

Figure 6: Fourier magnitudes in the (u, v)— plane, resulting from simulated observations 
with CTA layout B of binary stars with different diameters. Already small changes in the 
diameter of the secondary star by a few tens of /j,as show up clearly, also well visible in 
the image reconstructions of Figure [5] 
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images can be reconstructed on submilliarcsecond scales. A limitation that 
remains is the non-uniqueness between the image and its mirrored reflection. 

Figure [5] shows the results from such simulations of three binary stars, 
where the radius of one of them is varied. Already a change of the diameter 
of the secondary component by only a few tens of microarcseconds shows up 
clearly in the Fourier magnitude, and also the reconstructed images reproduce 
the stellar diameters and separations with quite satisfactory accuracy. 

The fidelity of the reconstructed image depends not only on 'obvious' 
factors such as the brightness of the target and the efficiency of the detectors 
but also on the position of the source in the sky, the geometric layout of 
the telescope array, and the latitude of the observatory. Aperture synthesis 
is achieved by the Earth's rotation carrying the star across the sky and - 
since the telescopes are fixed on the ground - the effective baselines, i.e., the 
separations between pairs of telescopes as seen along the line of sight, grad- 
ually change, filling in various portions of the (u, v)— plane. The geometry of 
the array and the celestial position of the source determine what projected 
baselines will be generated during the source's passage across the sky. For 
example, sources near the celestial poles do not move, and layout geometries 
with telescopes in repetitive patterns offer fewer unique baselines. 

7. 3. Stellar Diameters and Binary Separations 

The main purpose of the classical interferometer at Narrabri was to mea- 
sure angular diameters of stars, practical already with only two telescopes. 
It was also possible to study parameters such as binary separations by fitting 



models to the data [351 ] . With CTA, one will be able to perform such mea- 
surements in a much more accurate and model-independent manner, since 
such an array samples the Fourier plane in thousands of points as compared 
to the ~5 points for typical past measurements at Narrabri. For such stud- 
ies, utilizing some prior knowledge of the source (assuming it to be a binary 
star, for example), one can fit a model to the 'raw' data in the Fourier plane, 
without going into any (possibly algorithm-dependent) image reconstruction. 

7.4- Observations with subsets of the configurations 

Quite probably, not all CTA telescopes will be available for intensity in- 
terferometry at any one time. Even if most of the hardware installed for 
Cherenkov-light measurements could be utilized for stellar observations, its 
specific requirements with regard to signal handling (and possibly also auxil- 
iary optics and detectors), combined with finite resources, are likely to limit 
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Figure 7: Subsets of the candidate configurations (B, D, I from top to bottom). In the 
leftmost column one half of the telescopes of the superset configurations (Figure [2]) were 
selected in a pseudo-random fashion. In the middle column, one in four telescopes was 
selected, and in the rightmost column one in eight. 
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Figure 8: Simulated observations of the binary star in Figure |4] with subsets of configu- 
ration B (top row), D (middle row) and I (bottom row). The left column is for subsets 
containing half of the telescopes; the center for a quarter, and the rightmost for one eighth 
of the total (cf. Figure [7]). Best imaging is achieved with those arrays that have a balanced 
mix of different baselines. 
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the number of telescopes initially equipped for interferometry. Observations 
using only a subset of the telescopes may thus represent a realistic mode of 
operation, and we now consider the choice of subsets of telescopes. 

For each of the array layout configurations B, D and I (Figure [2]), three 
subsets were generated, as shown in Figure [7J The configurations shown in 
the leftmost column, designated B2, D2 and 12, were obtained by selecting 
half of the full set of telescopes in a semi-random manner, attempting to 
preserve the overall 'shape' of the array |54 |. In the middle column, one in 
four telescopes was retained (B4, D4, 14) and in the rightmost column only 
one telescope in eight was kept (B8, D8, and 18). 

Figure [HI shows the output from simulations of the binary star in Figure H] 
using these subsets. The magnitude of the star was now fixed to my = 5 and 
a long integration time was chosen in order to depress measurement noise 
and thus highlight sampling effects for the various configurations. 

It is obvious that more numerous telescopes, with a wider distribution of 
baselines, are better in terms of Fourier-plane sampling, and having only few 
telescopes restricts the results. Also, the optimal distribution of baselines 
depends on the actual size of the target. In this example, much of the infor- 
mation is at lower spatial frequencies, better sampled by the configuration B8 
which for those provides a denser sampling of the (u, v)— plane. Although D8 
in principle enables higher angular resolution, its sparse sampling with only 
a small variation of baseline lengths does not allow the Fourier magnitudes 
to be correctly estimated at low spatial frequencies, resulting in a somewhat 
blurred pattern. 

Note, however, that here it was not attempted to optimize the telescope 
selection for optimal sampling. 

8. The new stellar physics 

With optical imaging approaching resolutions of tens of microarcseconds 
(and with also a certain spectral resolution), we move into novel and pre- 
viously unexplored parameter domains. This requires attention not only to 
optimizing the instrumentation but also to a careful choice of targets to be 
selected which should be both astronomically interesting and realistic to ob- 
serve. With a foreseen brightness limit of perhaps my=6 or 7, and with 
sources of a sufficiently high brightness temperature, initial observing pro- 
grams have to focus on bright stars or stellar-like objects (26| . 
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Figure 9: Approximated angular diameters of the stars in the Bright Star Catalogue [471 ] . 
containing all stars in the sky with my < 6.5. Effective temperatures were estimated from 
the B — V color index listed there, using a polynomial fitted to the relation between that 
index and T e ff by Bessell et al. [8|. The temperatures were then used to calculate the 
angular diameters by approximating stars as blackbodies with uniform circular disks. 
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Among the about 9000 objects in the Bright Star Catalogue [47] 
2600 objects are both hotter than 9000 K and brighter than my= 7, among 
which the brightest and hottest should be those easiest to observe. A selection 
of some 35 stars brighter than my= 2 or hotter than T e ff = 25,000 K, and of 
special astrophysical interest were listed as candidates for early observations 



by Dravins et al. [28|, including the following categories: 



8.1. Rapidly rotating stars 

Rapid rotators are normally hot and young stars, of spectral types O, B, 
and A. Some are rotating so fast that the effective gravity in their equatorial 
regions becomes very small, easily enabling mass loss or the formation of 
circumstellar disks. Rapid rotation causes the star itself to become oblate, 
and induces gravity darkening: equatorial regions become dimmer, and polar 
ones brighter. 

A number of these have now been studied with amplitude interferome- 
ters. By measuring diameters at different position angles, the rotationally 
flattened shapes of the stellar disks are determined. For some stars, also their 
asymmetric brightness distribution across the surface is seen, confirming the 
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expected gravitational darkening and yielding the inclination of the rota- 
tional axes. Aperture synthesis has permitted the reconstruction of images 
using baselines up to some 300 m, correspo nding to resolutions of 0.5 mas in 



the near-infrared H-band around A 1.7 /xm [123 . 

Two stars illustrate different extremes: Achernar (a Eridani) is a highly 
deformed Be-star (V rot sin i = 250 km s -1 ; > 80 % of critical). Its disk is the 
flattest so far observed - the major/minor axis ratio being 1.56 (2.53 and 
1.62 mas, respectively); and this projected ratio is only a lower value - the 



actual one could be even more extreme [20j. Further, the rapid rotation of 



Achernar results in an outer envelope seemingly produced by a stellar wind 



emanating from the poles [57J, |58j. There is also a circumstellar disk with Ha- 
emission, possibly structured around a polar jet (Hi}. The presence of bright 
emission lines is especially interesting: since the S/N of an intensity inter- 
ferometer is independent of the spectral passband, studies in the continuum 
may be combined with observations centered at an emission line. 

Going to the other extreme, Vega (a Lyrae, AO V) has been employed as 
one of the standard northern hemisphere calibration stars for optical astron- 
omy, but its true spectrum has turned out to be quite complex. First, space 
observations revealed an excess flux in the far infrared, an apparent signa- 
ture of circumstellar dust. Later, optical amplitude interferometry showed 
an enormous (18-fold) drop in intensity at A 500 nm from stellar disk cen- 
ter to the limb, indicating that Vega is actually a very rapidly rotating star 
which just happens to be observed nearly pole-on. The true equatorial rota- 
tional velocity is estimated to 270 km s^ 1 ; while the projected one is only 22 
kms -1 (3I l95j]. The effective polar temperature is around 10,000 K, the equa- 
torial only 8,000 K. The difference in predicted ultraviolet flux between such 
a star seen equator-on, and pole-on, amounts to a factor five, obviously not 
a satisfactory state for a star that should have been a fundamental standard. 

8.2. Circumstellar disks 

Rapid rotation lowers the effective gravity near the stellar equator which 
enables centrifugally driven mass loss and the development of circumstellar 
structures. Be-stars make up a class of rapid rotators with dense equato- 
rial gas disks; the 'e' in 'Be' denotes the presence of emission in Ha and 
other lines. Observations indicate the coexistence of a dense equatorial disk 
with a variable stellar wind at higher latitudes, and the disks may evolve, 



develop and disappear over timescales of months or years [97J. The detailed 



mechanisms for producing such disks are not well understood, although the 
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Figure 10: Types of primary targets for kilomctcr-scale intensity interferometry. Top row: 
Stellar shapes and surfaces affected by rapid rotation - The measured sha pe of Achernar 
[20j j ; expected equatorial bulge and polar brightening of a very rapid rotator |l!3| ; deduced 
surface brightness of the rapidly rotating star Vega, seen pole-on (95[; possible donut- 
shape for a rapidly and differentially rotating star [75|- Middle row: Disks and winds 

- Modeled interferometric image of the circumstellar disk of the Be-star £ Tauri [l(|; a 
magnetic stellar wind compresses a circumstellar disk |97l|: simulation of how stronger 
magnetic fields distort wind outflow from hot stars |115| : the strongest stellar wind in 
a binary opens up cavities around the other star: the geometry around the Wolf-Rayet 
star 7 2 Velorum as deduced from interferometry [8l| . Bottom row: Stellar surroundings 

- Interferometric image of the giant star T Leporis surrounded by its molecular shell [72j ; 
an analogous image of the giant e Aurigae, while partially obscured by a circumstellar 
disk [61]; artists view of the interacting j3 Lyrae system with a gas stream, accretion disk, 
jet- like structures and scattering halo 43j; an adaptive-optics, high-resolution image of 
the mysterious object r\ Carinae, the most luminous star known in the Galaxy pSOl ]. 
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material in these decretion (mass-losing) disks seems to have been ejected 
from the star rather than accreted from an external medium. 



8. 3. Winds from hot stars 

The hottest and most massive stars (0-, B-, and Wolf-Rayet types) have 
strong and fast stellar winds that are radiatively driven by the strong pho- 
tospheric flux being absorbed or scattered in spectral lines formed in the 
denser wind regions. Not surprisingly, their complex time variability is not 
well understood. Stellar winds can create co-rotating structures in the cir- 
cumstellar flow in a way quite similar to what is observed in the solar wind. 
These structures have been suggested as responsible for discrete absorption 
components observed in ultraviolet P Cygni-type line spectra. 

Rapid stellar rotation causes higher temperatures near the stellar poles, 
and thus a greater radiative force is available there for locally accelerating 
the wind. In such the result may be a poleward deflection of wind 

streamlines, resulting in enhanced density and mass flux over the poles and a 
depletion around the equator (opposite to what one would perhaps 'naively' 
expect in a rapidly rotating star). Surface inhomogeneities such as cooler 
or hotter starspots cause the local radiation force over those to differ. This 
leads to locally faster or cooler stellar-wind streamers which may ultimately 
collide, forming co-rotating interaction regions. Further, effects of magnetic 
fields are likely to enter and - again analogous to the ca se of the solar wind 
- such may well channel the wind flow in complex ways jl!5| . 

8.4- Wolf-Rayet stars and their environments 

Being the closest and brightest Wolf-Rayet star, and residing in a binary 
jointly with a hot O-type star, r y 2 Velorum is an outstanding object for studies 
of circumstellar interactions. The dense Wolf-Rayet wind collides with the 
less dense but faster O-star wind, generating shocked collision zones, wind- 
blown cavities and eclipses of spectral lines emitted from a probably clumpy 



wind |8l|, |87|] . The bright emission lines enable studies in different passbands 



and already with the Narrabri interferometer, Hanbury Brown et al. [41 
could measure how the circumstellar emission region (seen in the C III-IV 
feature around A 465 nm) was much more extended than the continuum flux 
from the stellar photosphere, and seemed to fill much of the Roche lobe 
between the two components of the binary. 
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A few other binary Wolf-Rayet stars with colliding winds are bright 
enough to be realistic targets. One is WR 140 (my = 6.9, with bright emis- 
sion lines), where the hydrodynamic bow shock has been monitored with 
milliarcsecond resolution in the radio, using the Very Long Baseline Array 
(VLBA). This revealed how the bow-shaped shock front rotates as the orbit 
progresses during its 7.9 yr period [211 ] . 

8.5. Blue supergiants and related stars 

Luminous blue variables occupy positions in the Hertzsprung-Russell di- 
agram adjacent to those of Wolf-Rayet stars, and some of these objects are 
bright enough to be candidate targets, e.g., P Cyg (my =4.8). Luminous 
blue variables possess powerful stellar winds and are often believed to be 
the progenitors of nitrogen-rich WR-stars. Rigel (/3 Ori; B8 lab) is the 
closest blue supergiant (240 pc). It is a very dynamic object with variable 
absorption/emission lines and oscillations on many different timescales. Ac- 
tually, the properties of Rigel resemble those of the progenitor to supernova 
SN1987A. 

A most remarkable luminous blue variable is rj Carinae, the most luminous 
star known in the Galaxy. It is an extremely unstable and complex object 
which has undergone giant eruptions with huge mass ejections during past 
centuries. The mechanisms behind these eruptions are not understood but, 
like Rigel, rj Car may well be on the verge of exploding as a core-collapse 
supernova. Interferometric studies reveal asymmetries in the stellar winds 
with enhanced mass loss along the rotation axis, i.e., from the poles rather 
than from the equatorial regions, resulting from th e enhanc ed temperature 



at the poles that develops in rapidly rotating stars [1171 . Il 19 



8.6. Interacting binaries 

Numerous stars in close binaries undergo interactions involving mass flow, 
mass transfer and emission of highly energetic radiation, Indeed, many of 
the bright and variable X-ray sources in the sky belong to that category. 
However, to be a realistic target for intensity interferometry, they must also 
be optically bright, which typically means B-star systems. 

One well-studied interacting and eclipsing binary is Lyrae (Sheliak; 
my = 3.5). The system is seen close to edge-on and consists of a B7-type, 
Roche-lobe filling and mass-losing primary, and an early B-type mass-gaining 
secondary. This secondary appears to be embedded in a thick accretion disk 
with a bipolar jet seen in emission lines, causing a light-scattering halo above 
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its poles. The donor star was initially more massive than the secondary, but 
has now shrunk to about 3M , while the accreting star has reached some 
13 M . The continuing mass transfer causes the 13-day period to increase by 



about 20 seconds each year [43 



Using the CHARA interferometer with baselines up to 330 m, the /3 Lyr 



system has been resolved in the near- infrared H and K bands [122J. The 
images resolve both the donor star and the thick disk surrounding the mass 
gainer, 0.9 mas away. The donor star appears elongated, thus demonstrat- 
ing the photospheric tidal distortion due to Roche-lobe filling. Numerous 
other close binaries invite studies of mutual irradiation, tidal distortion, limb 
darkening, rotational distortion, gravity darkening, and oscillations. 

8. 7. Observing programs 

Promising targets for early intensity interferometry thus appear to be 
relatively bright and hot, single or binary 0-, B-, and WR-type stars with 
their various circumstellar emission-line structures. The expected diameters 
of their stellar disks are typically on the order of 0.2-0.5 mas and thus lie 
(somewhat) beyond what can be resolved with existing amplitude interfer- 
ometers. However, several of their outer envelopes or disks extend over a few 
mas and have already been resolved with existing facilities, thus confirm- 
ing their existence and providing hints on what types of features to expect 
when next pushing the resolution by another order of magnitude. Also, when 
observing at short wavelengths (and comparing to amplitude interferometer 
data in the infrared), one will normally observe to a different optical depth 
in the source, thus beginning to reveal also its three-dimensional structure. 

Also some classes of somewhat cooler objects are realistic targets. Some 
rapidly rotating A-type stars of temperatures around 10,000 K should be 
observable for their photospheric shapes (maybe even watching how the pro- 
jected shapes change with time, as the star moves in its binary orbit, or if 
the star precesses around its axis?). 

The exact amounts of observing time required for different targets are 
somewhat awkward to estimate since - in contrast to 'classical' observations, 
the achievable signal-to-noise ratio in intensity interferometry depends on 
several factors other than apparent magnitude (not least the source's own 
brightness temperature in either the continuum or in some spectral line). It 
is also a function of the (normally unknown) source structure: possible high- 
contrast features on the milliarcsecond level will produce more measurable 
Fourier power in the (u,v)— plane. Of course, the signal-to- noise improves 
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with higher detective quantum efficiency, better electronic time resolution, 
the number (and size) of telescopes used, and with the number of wavelength 
channels that are simultaneously handled; Eq.®. 

Various simulations (and an extrapolation from past experience with the 



Narrabri stellar intensity interferometer) have shown |28|, |71| that a few tens 
of hours of integration with a large array of the CTA type will faithfully re- 
produce the Fourier pattern in the (u, v)— plane for a star of visual magnitude 
m.y = 5, and T e ff — 10,000 K, when measured in one wavelength channel with 
a time resolution of 1 ns. From such data, a two-dimensional image clearly 
can be reconstructed. While significant progress on image reconstruction al- 



gorithms has recently been made [89l . |90| , it is not yet known how sensitive 



reconstructed images could be to various types of noise levels, either such due 
to limited source brightness, limited (u,v)— plane coverage, or instrumental 
systematics. If, in the event of noisy observations, the data would not per- 
mit full two-dimensional imaging, one could instead extract one- dimensional 
quantities such as the sizes of stellar surface structures, the amount of limb 
darkening, or the separation of binary components. 

9. Observing in practice 

In carrying out actual observations for intensity interferometry, various 
practical and technical issues may require attention, concerning aspects of 
the telescopes, detectors, data handling, and the scheduling of observations. 

9.1. Optical e-interferometry 

Electronic combination of signals from multiple telescopes is becoming 
common for long-baseline radio interferometry, where remote radio antennas 
are connected to a common signal-processing station via optical fiber links in 
so-called e-VLBI. This is feasible due to the relatively low radio frequencies 
(MHz-GHz); doing the same for a corresponding optical phase- resolved signal 
(THz-PHz) would not be possible but the much slower intensity-fluctuation 
signals (again MHz-GHz) are realistic to transmit, thus enabling an electronic 
connection of also optical telescopes. A number of authors have noted this 
potential of electronically combining multiple optical apertures, especially for 
observations at short optical wavelengths (perhaps using the multiple mirror 



segments of t he p rimary mirror in extremely large telescopes [23 



24|). Ofir 



& Ribak 9l|, |92j, |93( evaluate concepts for multidetector intensity interfer- 



ometers, and even space-based intensity interferometry has been proposed 
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52|, |60[, exploiting the possibility to combine signals off-line from each com- 
ponent telescope, thus relaxing the requirement for spacecraft orientation 
and orbital stability. 



9.2. Performance of Cherenkov telescopes 

The signals to be measured for intensity interferometry have much in 
common to those of atmospheric Cherenkov flashes: nanosecond time struc- 
ture and relatively short optical wavelengths. Most probably, the same types 
of very fast photon-counting detectors can be used, although the sources 
to be observed are much brighter, and the data handling has to allow for 
continuous integrations (rather than trigger-based acquisition of short data 
bursts). 

9.2.1. Image quality 

Even if the technique of intensity interferometry as such does not require 
good optical quality, and permits also rather coarse light collectors with 
point-spread functions of several arcminutes, issues arise from unsharp stellar 
images: in particular an increased contamination from the background light 
of the night sky. Although this light does not contribute any net intensity- 
correlation signal, it increases the photon-counting noise, especially when 
observing under moonlight conditions. 

While any reasonable optical quality should be adequate for intensity 
interferometry as such, the magnitude my of the faintest stars that can 
be studied will be influenced by the optical point spread function. Two 
extreme sky brightness situations can be: (a) dark observatory sky with 
~ 21.5 my/arcsec 2 ; (b) sky with full Moon; ~ 18my/arcsec 2 . The equivalent 
magnitudes from the sky background then result in my ~9.4 (a) and 5.9 (b) 
for a 5arcmin diameter field, and my ~ 12.9 (a) and 9.4 (b) for a 1 arcmin 
diameter field. 

A larger point spread function also takes in other sky events (meteors, 
distant flashes of lightning, etc.), and may preclude the use of small-sized 
semiconductor detectors of possibly higher quantum efficiency. 

9.2.2. Isochronous optics 

For Cherenkov light observations, a large field of view is desired. In 
most optical systems, the image quality deteriorates away from the optical 
axis, and to mitigate this, various optical solutions are used. Many current 
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telescopes have the layout introduced by Davies & Cotton [151 ]. whose pri- 
mary reflector forms a spherical structure centered on the focal point, giving 
smaller aberrations off the optical axis compared to a parabolic design. 

This has the consequence that the telescope optics become anisochronous, 
i.e., photons originally on the same wavefront, but striking different parts 
of the entrance aperture may not arrive to the focus at exactly the same 
time. As noted above, the signal-to-noise ratio improves with electronic 
bandwidth, i.e., the time resolution with which stellar intensity fluctuations 
can be measured. The time spread induced by anisochronous telescopes acts 
like 'instrumental broadening' in the time domain, filtering away the most 
rapid fluctuations. This probably is not a serious issue since the gamma- 
ray induced Cherenkov light flashes in air last only a few nanoseconds, and 
thus the performance of Cherenkov telescopes cannot be made much worse, 
lest they would lose sensitivity to their primary task. Still, since realistic 
electronics may reach resolutions on the order of 1 ns, it would be desirable 
that the error budget does not have components in excess of such a value. 

Among existing Cherenkov telescopes, this is satisfied by parabolic de- 
signs (e.g., MAGIC) but not by the Davies-Cotton concept (e.g., VERITAS or 
H.E.S.S.-I). For example, in the H.E.S.S.-I telesco pes t he photons are spread 
over At ~ 5 ns, with an rms width ~ 1.4 ns 0, @, 104 1. For large telescopes, 
the time spread would become unacceptably large if a Davies-Cotton design 
were chosen, and those therefore normally are parabolic (e.g., MAGIC on La 
Palma; H.E.S.S.-II in Namibia, and MACE in Ladakh, India). In principle, 
these then become isochronous - apart from minute (few hundred ps) effects 
caused by individual mirror facets being spherical rather than parabolic, or 
by the tesselated mirror facets being mounted somewhat staggered in depth. 

Also non-parabolic telescopes can be made effectively isochronous, if they 
have more than one optical element. The two- mirror Schwarzschild-Couder 
design is attractive for smaller telescopes, not least because its sma ller image 



scale permits smaller and less expensive focal-plane cameras [1181 ]. For on 



axis rays, this design in principle is isochronous, but the time spread increases 
to ~ 1 ns for angles a few degrees off center. Also, Schmidt-type telescope 
designs may satisfy high demands on isochronicity (even better than 10 ps 
on axis), while also being compact, offering a wide field of view, and having 
a narrow point-spread function j82|. However, to take full advantage of such 
performance would require corresponding accuracies in all other components 
of the error budget, including the signal handling, and the positioning of 
telescopes on millimeter scales. Such values also begin to approach the level 
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of natural fluc tuati ons in path-length differences induced by atmospheric 
turbulence [12|, Il21 . 



9.2.3. Focusing at 'infinity' 

The optical foci of Cherenkov telescopes are optimized to correspond to 
those heights in the atmosphere where most of the Cherenkov light originates, 
and the image of a distant star will then be slightly out of focus. For a focal 
length of / = 10 m, the focus shifts 1 cm between imaging at 10 km distance 
and at infinity, which for an f/1 telescope implies an additonal image spread 
of some cm. In order to decrease the stellar image and not to take in too much 
of the night-sky background, it is desirable (though not really mandatory) to 
refocus the telescope on stars at 'infinity'. On some (especially larger) tele- 
scopes, such a possibility may be available anyhow since some refocusing can 
be required in response to mechanical deflections when pointing in different 
elevations or as caused by nocturnal or seasonal temperature variations. In 
the absence of such a possibility, a refocusing could still be achieved by an 
optical lens placed directly in front of the photosensor. 

9.2.4- Placement of telescopes in an array 

The placement of telescopes in interferometers c an b e optimized for the 
best coverage of the (u,v)— plane 0, Hil 48, 59, 83, 112 ] . As the star grad- 



ually crosses the sky during the night, projected baselines between pairs of 
telescopes change, depending on the angle under which the star is observed. 
If the telescopes are placed in a regular geometric pattern, e.g., a repetitive 
square grid, the projected baselines are similar for many pairs of telescopes, 
and only a limited region of the (u,v)— plane is covered (on the other hand, 
redundant baselines result in better signal-to-noise for those particular ones). 
Since stars rise in the east, moving towards west, baselines between pairs 
of telescopes that are not oriented exactly east-west will trace out a wider 
variety of patterns. Because of such considerations, existing amplitude in- 
terferometers (both optical and radio) locate their component telescopes in 
some optimal manner (e.g., in a Y-shape, or in logarithmic spirals, unless 
constrained by local geography). 

As concerns specifically the CTA, its smaller telescopes will be so nu- 
merous that, for most practical purposes, their exact placement should not 
be critical for interferometry - a huge number of different baselines will be 
available anyway. However, the situation is different for the very few large 
telescopes. Avoiding placing them on a regular grid (such as a square) will 



31 



offer a variety of baseline lengths, give a better coverage of the (u,v)— plane, 
and permit better image reconstruction. 



9.2.5. Impact on observatory operations 

The impact of intensity interferometry on other Cherenkov array oper- 
ations should not be significant. An important aspect is that - while full 
moonlight may constrain observations of the feeble atmospheric Cherenkov 
light - measuring brighter stars is no problem for intensity interferometry, 
enabling efficient operations during both bright- and dark-Moon periods. 

Potential sources for interferometry are distributed over large parts of 
the sky and permit vigorous observing programs from both northern and 
southern sites. However, several among the hot and young stars belong to 
Gould's Belt, an approximately 30 million year old structure in the local 
Galaxy, sweeping across the constellations of Orion, Canis Major, Carina, 
Crux, Centaurus, and Scorpius, centered around right ascensions 5-7 hours, 
not far from the equator. Thus, many primary targets are suitable to ob- 
serve during northern-hemisphere winter or southern-hemisphere summer. 
We note that this part of the sky is far away from the many gamma-ray 
sources near the center of the Galaxy (which is at right ascension 18 hours). 

9.3. Detectors and cameras 

Typical Cherenkov telescopes have focal lengths on the order of 10 m, pro- 
viding a focal-plane image scale around 3 mm/arcmin. A typical point-spread 
function of 3 arcmin diameter thus corresponds to 1 cm. Detectors that are 
capable of photon counting with nanosecond time resolution include well- 
established vacuum-tube photomultipliers and large-size solid-state avalanche 
diode arrays that are under development. 

A Cherenkov telescope typically holds several hundred photomultiplier 
tubes acting as 'pixels' in its focal-plane camera. The detectors and their 
ensuing electronics are naturally optimized for the triggering on, and the 
recording of, faint and brief transients of Cherenkov light and might not 
be readily adaptable for hour-long continuous recordings of bright stellar 
sources. However, for intensity interferometry, only one pixel is required (at 
least in principle, although some provision for measuring the signal at zero 
baseline is required) and we note that in some telescopes (e.g., HEGRA [iij 
and MAGIC 74]|), the central camera pixel was specifically designed to be 



accessible for experiments without affecting any others. Such types of central 
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pixels could be usable to perform some experiments towards also intensity 
interferometry. 

However, even if a special pixel is accessible, it may not be possible to 
use it in its bare form. If observing a bright source in broadband white light 
with a large telescope, the photon-count rate may become too large to handle, 
even for reduced photomultiplier voltages. As discussed above, the signal-to- 
noise ratio in intensity interferometry is independent of the optical passband: 
the smaller photon flux in a narrow spectral segment is compensated by 
the increased temporal coherence of the more monochromatic light. This 
property can be exploited with some color filter to reduce the photon flux to 
a suitable level, or using a narrow-band filter tuned to some specific spectral 
feature of astrophysical significance. For such uses, there should be some 
provision for a mechanical mounting in front of the detector to hold some 
small optical element (s). A broader-band color filter could simply be placed 
immediately in front of a photomultiplier but a narrow-band filter could 
require additional arrangements. Most such filters are interferometric ones 
and need to be used in collimated (parallel) light in order to provide a more 
precise narrow passband. Since light reaching the Cherenkov camera is not 
collimated, some additional optics could be required, or else one might use 
narrow-band filters based on other optical principles, such as Christiansen 
filters 0. 

The further development and optimization of observational techniques is 
likely to involve experiments with other types of detectors, color filters, po- 
larizers or other optical components which could be awkward to mechanically 
and electronically (re)place in the regular Cherenkov camera. To minimize 
disturbances to the Cherenkov camera proper, it could then be preferable to 
place an independent detector unit on the outside of its camera shutter lid. 
Such constructions have already been made on existing Cherenkov telescopes, 
e.g., a 7-pixel unit on a H.E.S.S. telescope used a plane secondary mirror to 
put it into focus, and was used for experiments in very high time-resolution 
optical observations. Its central pixel recorded the light curve of the target, 
while a ring of six surrounding pixels monitored the sky back grou nd and 



acted as a veto system to reject atmospheric background events [17], [18|, |46 



For such devices, provision must also be made for electrical power supply 
and signal cables to/from the outside of these camera shutter lids. 
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9.4- Signal handling 

Current electronic units, used in various photon-counting experiments, 
have time resolutions approaching 1 ns, and the error budget should ide- 
ally not have components in excess of such a value (the signal-to-noise is 
proportional to the square root of the signal bandwidth). Telescopes may 
be separated by up to a kilometer or two, and the timing precision of the 
photon-pulse train from the detector to a central computing location should 
be assured to no worse than some nanosecond (for the timing of its leading 
pulse-edge; the pulse- width may be wider). Such performan ce ap pears to be 



achievable by analog signal transmission in optical fibers [98|, Il20j . Compared 
to metal cables, these have additional advantages of immunity to cross-talk 
and to electromagnetic interference, and also avoid the difficulty of maintain- 
ing a common ground and protection for the receiving electronics against (in 
some locations not uncommon) lightning strikes across the array. 

Another possibility is using clocks on satellite positioning systems such 
as GPS or Galileo, where absolute timing within some nanosecond has been 
achieved in astronomical instrumentation [84|. This enables to time-tag the 
photon stream for later off-line analysis, with an accuracy better than the 
anisochronicity of the mirror surfaces of the typical Davies-Cotton optical 
design. 

9.4-1- Correlators 

A critical element of an intensity interferometer is the correlator which 
provides the averaged product of the intensity fluctuations (AiiAI 2 ) to be 
normalized by the average intensities and (I2) (Eq. H|). The original in- 
terferometer at Narrabri used an analog correlator to multiply the photocur- 
rents from its phototubes, and significant efforts were made to shield the 
signal cables from outside disturbances. Current techniques, such as FPGA 
(Field Programmable Gate Arrays), permit to program electronic units into 
high-speed digital correlators with time resolutions of a few ns or better. 

Similar units are also commercially avail able for primary applications in 
light scattering against laboratory specimens |lOl| . Such intensity-correlation 
spectroscopy is the temporal analog to the [spatial] intensity interferometry, 
and was developed after its subsequent theoretical understanding. It was 
realized that high-speed photon correlation measurements were required and 
electronics initially developed in military laboratories were eventually com- 



mercialized, first by Malvern Instruments in the U.K. 96], and nowadays 
offered by various commercial companies Q. 
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An alternative approach (at least for limited photon-count rates) is to 
digitize and store all data, and then later perform the correlation analyses off- 
line. The data streams from multiple telescopes can then be cross-correlated 
using a software correlation algorithm, permitting the application of digi- 
tal filtering to eliminate possible interference noise from known sources, and 
also to compute other spatio-temporal parameters, such as higher-order cor- 
relations between three telescopes or more, which in principle may contain 
additional information. On the other hand, this requires a massive comput- 
ing effort and possible observational problems may not get detected while 
observations are in progress but only at some later time. Such a capability 
was foreseen in the design study for QuantEYE, a proposed very high-time 



resolution instrument for extremely large telescopes [24], |25j , and verified in 



the construction and operation of the AquEYE and IquEYE instruments, the 



latter used also at the European Southern Observatory in Chile [84j, |85j, |86 



9.5. Delay units 

Besides the correlator, another piece of electronics is required for real- 
time intensity interferometry, namely to implement a continuously variable 
time delay that compensates for the relative timing of the wavefront at the 
different telescopes, as the source moves across the sky (Eq. |7J). 

If such a delay unit is not used, the maximum correlation signal in a mul- 
tichannel digital correlator will appear not in the channel for zero time delay 
between any pair of telescopes, but rather at that channel which corresponds 
to a delay equal to the light-time difference between telescopes along the 
line of sight towards the source. This arrangement is feasible already with 
existing digital correlators since these can be programmed to measure the 
correlation at full time resolution at time coordinates away from zero. 

Such arrangements, however, are not required in the case of off-line data 
analysis, where the delays can be introduced by software afterwards. 

10. Experimental work 

As preparatory steps towards realizing full-scale stellar intensity interfer- 
ometry, different laboratory and field experiments have been carried o ut a t 



various institutes, in particular at the StarBase facility in Utah |28|, [68|, 1109 
Also, in a first full-scale test with a Cherenkov telescope array, pairs of the 
12 m telescopes of the VERITAS array in Arizona were used to observe a 
number of stars, with pairs of its telescopes interconnected through digital 
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correlators 27]. For these observations, starlight was detected by a photon- 
counting photomultiplier in the central pixel of the regular Cherenkov- light 
camera, the outgoing photon pulses were digitized using a discriminator, then 
pulse-shaped and transmitted from each telescope via an optical cable to the 
control building where they entered a real-time digital cross correlator, com- 
puting the cross correlation function for various time delays. Continuous 
count rates up to some 30 MHz were handled, limited by the digitization 
and signal-shaping electronics. While these experiments were not intended 
to measure astrophysical quantities but to gain experience in operating with 
a full-scale observatory, they confirmed that no fundamental problems seem 
to exist in carrying out such operations. 

11. Further possibilities 

The availability of very large light-collecting areas, distributed over an 
extended array enables further classes of optical observations, not feasible 
with ordinary instrumentation. 

11.1. Higher- order spatio-temporal correlations 



The quantum theory of optical coherence |33|, |76| describes how one can 
define correlations between arbitrarily many spatial and/or temporal coordi- 
nates in the volume of light ('photon gas') being received from a source. The 
spatial intensity interferometer is only one special case of such more gen- 
eral spatio-temporal correlations, in that it measures the cross correlation 
between the intensities at two spatial locations, at one instant in time. 

However, using telescope arrays, and given that their photon detectors 
provide data streams which can be analyzed at will, one can construct, 
e.g., third-order intensity correlations, for systems of three telescopes: 
(I(ri, ti)I(r 2 , t 2 )I(r 3 , t 3 )), where the temporal coordinates do not necessarily 
have to be equal. In principle, such and other higher-order spatio-temporal 
correlations in light may carry additional information about the source from 
where the light has been emitted and thus - at least in principle - is of rele- 
vance for astronomy where information about the source has to be extracted 



from more or less subtle properties of its radiation received [53J, |91 . 

Although, in the recording of higher-order correlations, also the relative 
noise level increases (possibly demanding very large telescopes for certain 
measurements (221), all sorts of higher-order correlations can in principle be 
obtained without any additional observational effort if the digital signals from 
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each telescope are available for further manipulation in either hard- or soft- 
ware. For example, one could calculate correlations among also all possible 
triplets and quadruplets of telescopes, possibl y enabling a more rob ust full 
reconstruction of the source image 0, 0, HH Eq3, H, H, [l24 . 



11.2. High-time-resolution astrophysics 

Further uses of CTA can be envisioned in searching for extremely rapid 
(micro- or even nanosecond) optical variability, such as suspected from pul- 
sars or other compact sources. In a few radio pulsars, nanosecond pulse 
structure has been observed, representing the currently most rapid fluctua- 
tions seen in astronomical sources but so far lacking any credible explana- 
tion (perhaps nonlinear plasma turbulence, stimulated Compton scattering, 
or angular beaming due to relativistic source motion?), and the re is some 
evidence that corresponding events may exist also in the optical jl07 |. De- 



spite the point-spread function extending over arcminutes (and thus giving 
a considerable contribution from the night-sky background), the very large 
collecting area distributed over several independent telescopes would make 
such searches more sensitive than with any existing large telescope. Although 
the number of photons per second in the background skylight may be very 
significant, their number per microsecond is still very modest and the sensi- 
tivity to detecting the very shortest fluctuations becomes mainly a function 
of light-collecting area [66J . 

The sensitivity of Cherenkov telescopes to detecting such very brief op- 
tical flashes was analyzed by Deil et al. [l8j]. Comparing a H.E.S.S. 100m 2 
Cherenkov telescope with a sky-background-free optical telescope for very 
high-time-resolution photometry, they found that for flashes shorter than 
some 100 ns, a large Cherenkov telescope outperforms today's largest astro- 
nomical telescopes, at least under dark-Moon sky conditions. Optical ob- 
servations of the Crab pulsar have been made with Cherenkov telescopes of 



H.E.S.S. |46| and MAGIC |74|. Although this optical pulse is only some 10~ 6 
of the background (the sky plus the Crab nebula surrounding the pulsar), 
and integrations over many hundreds of pulsar periods are required before 
a sensible signal appears, statistical information on also the very fine time 
structure can be retrieved within very reasonable integration times. 
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12. Outlook 



Interferometry for the attempted measurement of stellar diameter s ap 



pears to have been first carried out in the 1870's by Edouard Stephan 110 



following a suggestion by Hippolyte Fizeau. A two-aperture mask was placed 
over a 80 cm reflector at Marseille Observatory, but it was soon realized that 
stars could not be resolved over this short baseline. In the 1920's, Albert 
Michelson & Francis Pease 79J operated a 6-meter interferometer mounted 
on the 100-inch Hooker telescope on Mt. Wilson, and succeeded in measuring 
diameters of a few giant stars, while their later 15-meter instrument proved 



mechanically too unstable for practical use [42 



The demanding requirement to maintain stable optical path differences 
during observations within a fraction of an optical wavelength caused the 



suc- 



technique to lay dormant for half a century, until Antoine Labeyrie [62 
ceeded in measuring interference fringes between two separated telescopes. 
This success triggered the construction of a whole generation of optical am- 
plitude interferometers and is also said to have been the specific reason why 
the plans to build a successor to the original Narrabri intensity interferometer 
(designed around that very time) were not realized, and (as far as astronomy 
is concerned), the technique has now been dormant for decades. 

However, the progress in instrumentation and computing technology since 
the days of the Narrabri interferometer has been extraordinary. High-speed 
photon-counting detectors and hardware correlators are commercially avail- 
able, and new mathematical algorithms allow for image reconstruction. The 
most valuable components - large light collectors - are being realized in the 
form of air Cherenkov telescopes. All of this has sparked a renewed inter- 
est in astronomical intensity interferometry, and a first workshop (since very 
many years) on this topic was held not long ago [67|. Thus, long after the 
pioneering experiments by Hanbury Brown and Twiss, the technological de- 
velopments carry the promise of achieving a basic but difficult goal: to finally 
be able to view our neighboring stars not only as mere unresolved points of 
light but as the extended and most probably very fascinating objects that 
they really are. 
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